Mapping the braiding properties of the Moore-Read state 
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In this paper we explore the braiding properties of the Moore-Read fractional Hall sequence, 
which amounts to computing the adiabatic evolution of the Hall liquid when the anyons are moved 
along various trajectories. In this work, the anyons are pinned to precise spatial configurations by 
using specific external potentials. Such external potentials break the translational symmetry and it 
appears that one will be forced to simulate the braidings on the entire many-body Hilbert space, an 
absolutely prohibitive scenario. We demonstrate how to overcome this difficulty and obtain the exact 
braidings for fairly large Hall systems. For this, we show that the incompressible state of a general 
(A;, m) fractional Hall sequence can be viewed as the unique zero mode of a specific Hamiltonian 
^{k,m) ^ whose form is explicitly derived by using k-particles creation operators. The compressible 
Hall states corresponding to nxk anyons fixed at wi,. . . ,Wnk are shown to be the zero modes of a 
pinning Hamiltonian Hwi^.^w^^^ which is also explicitly derived. The zero modes of Hwil^.\w^f. are 
shown to be contained in the space of the zero modes of Therefore, the computation of the 

braidings can be done entirely within this space, which we map out for a number of Hall systems. 
Using this efficient computational method, we study various properties of the Moore-Read states. 
In particular, we give direct confirmation of their topological and non-abelian properties that were 
previously implied from the underlying Conformal Field Theory (CFT) structure of the Moore-Read 
state. 

PACS numbers: 71.10.-w, 71.15.-m 



I. INTRODUCTION 

The recently discovered Fractional Hall states at filling 
factors u = 5/2,— and z/=12/5,^ are thought to contain 
anyons obeying Non-Abelian statistics. It was arguecP 
that these states are part of the Moore-Read sequence^ 
and of the Read-Rezayi sequence at /c=3,^ respectively. 
Both states received sustained attention from the the- 
oretical and experimental condensed matter community 
because they can provide the key to the physical real- 
ization of a Topological Quantum Computer (TQC). For 
this reason, there is a concentrated effort on both theo- 
retical and experimental fronts for finding direct confir- 
mations of the Non-Abelian and Topological properties 
of these Fractional Hall states .'^'Hl 

In spite of an overwhelming indirect evidence, the non- 
Abelian statistics for the Moore-Read sequence was not 
considered resolved^^ until the relatively recent numer- 
ical confirmation of Ref. Ql^ where Monte-Carlo tech- 
niques were used to compute a small number of adiabatic 
mononodromies. The present paper reports the results of 
an exact diagonalization study, which allows us to take a 
direct and unprecedented look into the non-Abelian and 
topological properties of the Moore-Read sequence. We 
have results for systems containing 2 and 4 anyons. A 
model Hamiltonian (reported here for the first time) is 
used to pin down the anyons at specified locations and 
to move them adiabatically along different paths. Among 
our results, the reader will find the following: 

- a map of the Abelian and non-Abelian adiabatic cur- 
vature experienced by an itinerant anyon when the rest 
of the anyons are kept at fixed locations. 

- a map of a newly introduced Twist Density, which 



measures the twist of the zero modes space during the 
adiabatic braidings. 

- a direct proof that the adiabatic curvature is strongly 
localized near the fixed anyons, thus confirming the topo- 
logical properties of adiabatic braids. 

- a direct proof that the monodromy corresponding to a 
braiding in which one anyon loops around another anyon 
is in perfect agreement with the CFT prediction. 

- a direct proof that the zero modes space splits ac- 
cording to the fusion rules of the underling Conformal 
Field Theory structure of the Moore-Read state. 

- a direct proof that the conformal blocks can be dis- 
tinguished by bringing two anyons together and by mea- 
suring the electron density for the fused anyons. 

Let us briefly mention the relevance of our study for 
for TQC. The logic gates in TQC algorithms can be gen- 
erated by sequences of braids. The specific sequences 
of braids can be computed using the Solovay-Kitaev 
algorithm, once a unitary representation of the anyons' 
braid group is provided. For the Fibonacci anyons, for ex- 
ample, the braiding sequences implementing the elemen- 
tary logic gates were explicitely calculated in Refs.'Wand 
[TTl The unitary representations of the braid group can be 
classified by various methods. Particularly, the Ref. [18] 
describes a general method for generating unitary rep- 
resentations of the braid group using quantum groups. 
These representations were connected to the braidings 
of the anyons in Fractional Hall Liquids via the argu- 
ment that braiding can be computed from an effective 
picture in which the anyons carry an internal quantum 
symmetry. This symmetry can be derived from the un- 
derlying CFT structure of the Fractional Hall sequence. 
For example, the braiding of the Fibonacci anyons in 
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the sl{2)^ WZW theory can be exphcitly generated us- 
ing the Uq{sl{2)) quantum group at the root of unity 
^_g27r*/(/c+2) 18 ^g^g g^igQ arguecP^ that the braiding 
of the Read-Rezayi parafermions can be exphcitly gen- 
erated using the Uq^{sl{2)) (g) CZ4k,q2 quantum group at 
the roots of unity ^-^=6^^*/^^+^) and q2=e~'^^^'^^ . For 
the Moore-Read sequence (/c=2), the braiding reduces to 
the spinor representation of the mirror reflections in the 
group S0(2n), as discovered long time ago by Nayak and 
Wilczek.19 

In all these algebraic derivations the anyons are as- 
sumed infinitely far apart from each other. In practice 
this will not be the case, and this is the main reason why 
the algebraic results need a direct confirmation showing 
that the braiding properties remain unchanged when fi- 
nite anyon densities are considered. So far, our numerical 
results confirm that this is the case. 



These incompressible states occur at electron densities 
given in Eq. ([T]). 

Let us elaborate on the above conditions. From the 
definition given in Eq. [2j we have 

[vT'^^'^'^'^^]{zS,...,ZN) = Jd^ZlJd^Z2 

x(zi* +4)^(^1* -4)™'e-j(l^^l'+l^=^l') (4) 
x*('='™)(zi,^2,...,^iv)- 

The general structure of ^l/e^'™) is 

*('='™)(^l,...,^yv) = *fc(zi,...,Zjv) 

X n (^i-z,-re-i2:r=,N.p, (s) 

i<j<N 



II. FRACTIONAL HALL SEQUENCES: 
GENERAL CONSIDERATIONS. 



The Hall sequences can be label by two integers: 

A:=l,2,. . . , and m=0,l, The filling factor and the 

conformal charge for each sequence are given by: 



k 



A:m + 2' 



1 



2{k-l) 
k^2 



(1) 



The /c=l, u=l/{m + 2) sequence corresponds to the 
Laughlin state,^^ which has fractionally charged anyons 
carrying integer fiux. Their braiding is Abelian. The 
k>2 sequences have fractionally charged anyons carrying 
fractional fiux. Their braiding is believed to be Non- 
Abelian. The case k=2 corresponds to the Moore-Read 
sequence;^ k=3 and greater correspond to the Read- 
Rezayi sequences.^ Excepting /c=l, 2 and 4, all sequences 
support universal quantum computation. 

We define the /c-particle creation operators, t?^'^^^ 
generating (from the vacuum) /c-particles in a Laughlin 
state (we set the magnetic length to 1): 



X n 

i<j<k 



= (zi + ...+Zfe)^ 



(2) 



where z=x-\-iy is the complex representation of the posi- 
tion r={x^y) in the 2-dimensional plane, and M and m 
are integers larger or equal to zero. The Fractional Hall 
sequences can be described as follows. 

Incompressible states. We claim that the incompress- 
ible Hall state for an arbitrary (k^m)^ originally defined 
in terms of certain correlators of the Z/^ parafermion con- 
formal algebra, is the highest electron density state sat- 
isfying 



^(2,m )|^(fc,m)^ = 0, V M > and < 

^(fc+l,m)|^(fe,m)^ VM>0. 



(3) 



where ^j^{zi,...,zn) is the correlation function of N 
number of parafermion fields (we include in the 
factors that make the correlation function non-singular). 
We shall see that the first condition of Eq. (|3| relates 
to the part of the wavefunction contained in the second 
row of the above equation. Indeed, the right hand side 
of Eq. Q reduces to 

2<i<j 



x(^l-^2)" U{^^-^ln^^-^2r' 
i=3 



(6) 



X§/e(zi, ...,Zn) 

The multiple integral projects the part of the integrand 
that is invariant to 4-dimensional rotations in the (ri,r2) 
sub-space and we argue that this part is zero. We focus to 
rotations around the axis connecting ri+r2 to the origin. 
The second row of the above equation is invariant to such 
rotations but the invariant part of the remaining rows is 
identically zero. Indeed, the third and fourth rows are 
analytic functions in all complex variables, hence their 
product admit an expansion like Ca{zi — 2^2)", where 
Ca may depend on zi-\- Z2, zs^ ... ^zn- The invariant part 
of such expansion is equal to cq, but cq is identically zero 
when m'<m. This proves the first condition of Eq. [s] 
We consider now the second condition, which says that 



x(^r + ...,+4+i)^ n {z*-z*re 



fc+1 



i<j<A;+l 



X*(fc.™)(2i,...,2^). 



(7) 



is identically zero. The right hand side of Eq. ^ is equal 
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to: 



n {zi-Zj)'^e ' Jd^Zi...Jd^Zk+i 

k+2<i<j 



fc+i 



i<j</c+l 

x*fe(^i,...,zAr) n n (^i-^j-r- 

i=l j>/c+l 

(8) 

The third row is an analytic function in ah complex ar- 
guments, hence admits an expansion such as: 



OL i<j<k-\-l 



(9) 



where the coefficients may depend on 2:1 + ... + Zk-\-i^ 
Zk-\-2, . . . ,ZAr. The index a refers to the collection of in- 
dices aij appearing inside the product. The only nonzero 
contribution to the integral of Eq. [8] comes from the term 
a = 0, as one can see by switching to the coordinates 
s = (zi + . . . + Zk-\-i)/{k -\- 1)^ Zi = Zi — s. But there is 
no such term because ^/^(zi, . . . , zn) cancels identically 
when k-\-l particles come to the same point in space. ^ 

We reformulate Eq. [3] in terms of the creation oper- 
ators r]^^'^\wy generating /c-particle Laughlin clusters 
centered at different w's: 



(ri,...,rfc|7/(^'-)Ht|0) = 

n {z,-z,re 



(10) 



The incompressible states are then uniquely defined as 
the highest electron density states satisfying 



(11) 



Compressible states. The lower electron density states 
satisfy the same conditions but the uniqueness is, of 



course, lost. A lower density state "^w^ 



with n X k 



anyons located at {w}=wi^. . . ^Wnk has the general form 

i<j<N 

where . . . , zn) is a certain correlator.!^ This cor- 

relator is a multi- valued function in w^s. Each branch 
defines a different state. All these states have the anyons 
fixed at the same locations {w}. The degeneracy of the 
states with fixed anyon positions is discussed in details in 
Ref.fTSl for arbitrary k. The correlators ^|^|(zi,...,ZAr) 
have the property of vanishing whenever k particles meet 
at any of the w^s^^ 



We now show that these states can be defined by an 
additional condition to Eq. (|3|, namely: 



(13) 



(14) 



Indeed, we have 

[r?('='™)(«;^)*St;^)](^fe+i,...,^^) = 

k 

— - V \z-—w P 

jd^z,...jd^zk n (z*-z*re • 

^ N 2 

x^{^y{zi,...,ZN) n {Zi-Zj)'^e ' . 

i<j<N 

The right hand side can be written as: 

n (z,-z,-)'"e"'«^'='"'' Jd^zi...fd^Zk 

k<i<j 

X n Izi-Zj^e ^=1 

i<j<k 

k 

X§|^|(zi,...,ZAr) n Ui^i-^j)"^- 

i—lj>k 



The integral can be show to be identically zero by using 
the clustering properties of the correlator. 

Of course, we also need to prove that Eqs. ([3| and ([l3| 
are not satisfied by states other than the quasi-hole states 
of Eq. (12). For k=2 we have verified numerically (in the 
sphere geometry) that the number of linearly indepen- 
dent states satisfying these conditions is precisely equal 
to 2^~^, as it should. A similar study was conducted 
for the k = 3 Hall sequence. 



(15) 



III. THE PINNING HAMILTONIAN 

Based on the description given in the previous Sec- 
tion, we can generate simple model Hamiltonians for each 
Fractional Hall sequence. Such a task was previously un- 
dertaken in Ref. 25, 

We claim that all incompressible fractional Hall states 
corresponding to given k and m are zero energy states 
for the following Hamiltonian: 



(16) 



where A's must be all positive. For k=2 we have also 
verified that the opposite is true, namely that the zero 
modes space of the above Hamiltonian contains all the 
fractional Hall states and that its dimension coincide with 
the theoretical value derived in Ref. 23, A similar study 
exists for k=3 sequence.^ 
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FIG. 1: (Color online) a) and a') The particle density and b) 
and b') the pair amplitude corresponding to the unique zero 
mode of the pinning Hamiltonian H{wi,W2)[no- The first 
row corresponds to the odd or S=l sector (A^=15, A^0=29) 
and the second row corresponds to the even or S=0 sector 
(A/'=16, A/'(^=31). The position of the probes is indicated by 
red dots. 



states of the pinning Hamiltonian, for arbitrary anyon 
configuration {w} = wi^. . . ^Wnk^ can be computed in two 
steps without involving any approximation. Here is how: 

1. Construct the null space Ho of the Hamiltonian 



2. Restrict the pinning Hamiltonian: 

nk 

H{{w}) = Y,V^^^'^\w^yr]^^^^\w^) (18) 
/i=i 

to Ho and construct its null space l-Lo{{w}), which 
contains all the Hall states with anyons pinned at 

As we shall exemplify later, the dimension of the many- 
body Hilbert space increases extremely fast with the 
number of electrons. For this reason, any attempt of 
direct diagonalization of the full pinning Hamiltonian on 
the full Hilbert space is futile. Now, the difference be- 
tween the pinning Hamiltonian and is that the lat- 
est is translational invariant. In the sphere geometry,!^ 
this means that commutes with the total angular 
momentum L. Thus, to generate the null space Ho in 
the sphere geometry, we need to search only for the zero 
modes of highest weight, i.e. the ones in the null space 
of L+ operator. Once this step is completed, we can gen- 
erated the full Ho by successively applying L~ operator 
on the highest weight zero modes. This program was 
numerically implemented using standard techniques. 



IV. RESULTS FOR THE MOORE-READ 
SEQUENCE 



We also claim that a lower density state with nxk 
anyons present at Wi^. . . ^Wnk is a zero energy state of 
the following pinning Hamiltonian: 

hL^^::Iu,„, = xj£w ^(fc+i.™)(«;)tr,(fe+i.™)(ti;) 

+ J2 A^./d^ti; ,;(2,™')(y;)t^(2,™')(y;) 

m' <m (l'^) 

nk 

+ E A;. 7?(^'-)(^^)t7y(^'-)(^^), 

where, again, all A's are considered positive. It is useful 
to think of the pinning Hamiltonian H{{w}) as simulat- 
ing a set of nxk external probes that have been placed 
at the locations tL'i, . . . , Wnk- Again, for /c=2 and small 
number of units of flux added to the fundamental value 
of the magnetic flux, we have verified numerically (on 
the sphere) that the dimension of the zero energy modes 
space of the above Hamiltonian is precisely equal to 2^~^. 
A similar study exists for the /c=3 sequence. 

Given the particular form of the pinning Hamiltonian, 
namely, the fact that each term is positive definite, we 
have the following crucial observation: the zero energy 



Here we show a set of results for the Moore-Read se- 
quence /c=2 and m=l in the first Landau level. The 
calculations were performed with the sphere geometry, 
where we can work with a finite number of electrons. 
The number of available orbitals in the first Landau level 
is A^orb=A^(/)+l, where A^^ is the magnetic flux (in the 
quantum units of flux) passing through the surface of 
the sphere. We have to separately discuss the case of 
odd and even number of electrons. The incompressible 
state exists only for even number N of electrons, and oc- 
curs when the number of flux is relation N^=2N—3 with 
the number N of electrons. For this case, dim(Ho)=l- 

A. One pair of anyons. 

When the number of flux becomes A/'0=A'^+1, a pair 
of anyons is generated and 

dim(?^o) = lm2] + l){[N/2] + 2), (19) 

where the square brackets indicates the integer part.l^ 
We have computed and saved the Ho spaces for 
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FIG. 2: (Color online) The spectrum of H(wi, W2) [no func- 
tion of distance between the probes, which are moved along 
the meridians {0, = 0) and {0, = tt), with increasing from 
to 7r/2. The left /right column corresponds to odd/even 
number of electrons. Starting from the top, the left panels 
correspond to Ar/Ar^=9/16, 11/20, 13/24, 15/28 and the right 
panels to iV/iV^=10/18, 12/22, 14/26, 16/30. 



N/N^=4./6, 5/8, 6/10, 7/12, 8/14, 9/16, 10/18, 11/20, 
12/22, 13/24, 14/26, 15/28 and 16/30. The largest sys- 
tem in this sequence, the 16 electrons on 31 orbitals, has 
dim(Ho)=45. 

If we fix the positions of the two quasi-holes at arbi- 
trary locations Wi and by adding the pinning poten- 



tial described in Eq. 18, all the states in Ho are pushed 
up in energy, except one state whose energy remains ex- 
actly zero. This is precisely the state ^^/^'^^ discussed 
above. To exemplify, we consider the largest systems we 
computed so far, namely, the system with 15 electrons on 
29 orbitals (odd number of electrons) and the system with 
16 electrons on 30 orbitals (even number of electrons). 



In this cases, the total many-body Hilbert spaces have 
staggering dimensions of 77,558,760 and 300,540,195, re- 
spectively. In this extremely large Hilbert spaces, we 
find a number of zero modes for i^^^'^) equal to 36 in 
the first case and 45 in the second case (in total agree- 
ment with Eq. ([l9|). If we fix wi and W2 and diagonal- 
ize H{wi^ W2) [uoi we find one zero mode for both cases. 
Finding these zero modes would have been impossible 
without taking full advantage of the translational sym- 
metry at the first step, when Ho was resolved. 

To visualize a state, we compute the corresponding 
particle density and pair amplitude as functions of posi- 
tion on the sphere. The latest is given by the expectation 
value of iri^'^''^\wy iri^'^''^\w) on the zero mode. A plot 
of these quantities for the zero modes discussed above, 
is shown in Fig. [l] The positions of the probes were 
chosen as (6>=7r/2,^=0) and (6>=7r/2,^ = tt), so that we 
have maximum possible separation between the trapped 
anyons. Since the anyons are far apart, there is no vis- 
ible difference between even and odd cases. Refering to 
the Bratelli diagram for the Moore-Read sequence, the 
even and odd number of electrons correspond to the q- 
spin S'=0 and 5=1 sectors, respectively. Thus, we can see 
that there is no difference in the local properties of the 
wavefunctions belonging to different conformal blocks, 
which is precisely what one should see for a topologi- 
cal degeneracy. As we shall see, things look completely 
different when we bring the anyons close to each other. 
Other things to notice about Fig. [T] are the fact that the 
density is finite while the pair amplitude is exactly zero 
at the probe locations and the fact that the two anyons 
appear to be totally separated. 

Let us take a few lines here and explain our plots. 
Quantities that depend on the position on the sphere will 
be shown as surface plots, with the quantity of interest 
on the z axis. The cartesian coordinates x and y describe 
points of the sphere. If 9 and (j) are the usual angles on 
the sphere, then the relation between (6>,0) and (x, y) is 
given by 0=^ x'^ + y'^ and (/)=arctan(?//x). 

Next, let us take a look at the energy spectrum of 
H{wi^W2)[nQ cis function of probe separation, d{9) = 
^^N(f) sin I , while gradually increasing the number of elec- 
trons from 9 to 16. For each size, the probes were moved 
along the meridians (6>,^=0) and (^,^=7r), with increas- 
ing from to 7r/2. The strength of the probe potential 
was fixed at A=l. The results are shown in Fig.|2j where 
each panel displays a number of bands (equal to 
dim (Ho)), representing the flow of the eigenvalues with 
the distance d. There is one and only one eigenvalue 
that remains strictly zero (within a numerical error that 
is less than 10~^^!). The energy gap separating this zero 
mode from the rest of the spectrum goes to zero as the 
probes come closer to each other and converge to a well 
defined limit when the probes are moved far apart from 
each other. Compared to the other level spacings in the 
graph, the energy gap appears large. Another thing to 
notice is that there is a difference in the eigenvalues flow 
patterns when the plots for odd and even number of elec- 
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FIG. 3: (Color online) The particle density for fused probes 
as function of distance from the fusing point. Left /right 
panel refers to odd/even number of electrons. On the left, 
the different curves correspond to N/N^^9/IQ, 11/20, 13/24, 
15/28 and, on the right, the different curves correspond to 
iV/iV^=10/18, 12/22, 14/26, 16/30. 



trons are compared. It is also interesting to remark that 
the flow patterns remains almost unchanged as the size 
of the system is increased. 

Fussing the anyons. As one can clearly see in Fig. |2j 
when the anyons are fused, the lowest energy level be- 
comes degenerate. The electron density is not well de- 
fined, but we can still study its limit as the distance be- 
tween the anyons go to zero. In this limit, the electron 
density becomes radially symmetric and we can plot the 
density as a function of the distance from the position 
of the fused anyons. The results are shown in Fig. [3] for 
different sizes of the Hall system. The graphs reveal that 
the particle density is different for odd and even num- 
ber of electrons. This means that one should be able 
to tell when a pair of anyons is in the S=0 or S=l sec- 
tor by simply fusing the anyons together and measuring 
the electron density. Fig. |3] also reveals that the electron 
density is rapidly converging with the size of the system. 



B. Two pairs of anyons 

When the number of flux becomes A/'^ 
pairs of anyons are created anc^ 



=M+2, two 



dim(Ho) 



192 \ 



{N- 
{N- 



1) {N- 

2) (iV- 



■3)(A^4 
-4)2(iV- 



5)(iV- 
-6), 



■7) 



(20) 

where the first row refers to odd and the second row to 
even number of electrons. We have computed and saved 
the Ho spaces for N/N^=A/7, 5/9, 6/11, 7/13, 8/15, 
9/17, 10/19, 11/21, 12/23, 13/25, 14/27, 15/29 (16/31 
resisted to us so far). The largest system in this sequence, 
the 15 electrons on 30 orbitals, has dim(Ho)=660. 

It is instructive to consider an example. We pick 
the largest system we could compute so far, the A^=15 
and A/orb=30 case, where the total Hilbert space has 
a dimension of 155,117,520. In this extremely large 
Hilbert space, we find 660 zero modes for i^^^'^) (as 
we should). If we fix il'i,. . . , and compute the zero 
modes of H{wi^ . . . , w^) [^^^ we find two zero modes (as 




(a') 



(b') 




FIG. 4: (Color online) Visualization of the two zero modes 
for A/A0= 15/29. (a) and (a') shows the electron density 
and (b) and (b') shows the pair amplitude of the two zero 
modes. The zero modes were obtained by splitting the two 
dimensional space. The probes were fixed in a tetrahedral 
configuration. 



we should) . Fig. [4] shows the density and the pair am- 
plitudes corresponding to an orthogonal splitting of the 
zero modes space. The probes were fixed in a tetrahedral 
configuration for maximum separation. The orthogonal 
splitting of the zero modes space was generated by adding 
an infinitesimal, degeneracy lifting. Coulomb interaction 
to H{wi, . . . , W4). The position of the probes is indicated 
by red dots. One can see that, at the probe locations, 
the electron density is finite while the pair amplitude is 
exactly zero. Also, judging by these plots, one could say 
that the quasi-holes are separated. We should also point 
out that plots for the two zero modes look very similar, 
which is again the signature of the topological degener- 
acy. 

We continue our analysis of the case A/'0=A/'^+2 by 
plotting the spectrum of H{wi^ . . . , W4) for different 
probe locations. For this, we moved the probes continu- 
ously along the paths shown in the inset of Fig. [5j which 
can be described by: wi{0) = {0^(j)=0)^ W2{0)={0^(j)=7r\ 
^3(6>)=(7r-6>,0=7r/2) and ^4(<9)=(7r-6>,0=-7r/2). Fig. ^ 
plots the eigenvalues of H{wi, . . . , W4) [q-i^ as functions of 
6>, for N=1A and A'0=27 (the graph becomes extremely 
busy if larger systems are used). For 9 different from 
or TT, one can see in Fig. [5] the doubly degenerate zero 
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FIG. 5: (Color online) The spectrum of H{wi , . . . , wa) [no as a 
function of the probes' position. The probes were moved along 
the for meridians: (^,0 = 0), (^,0 = tt), (tt — ^, = vr/2), 
(tt — 0,(j) = — 7r/2) (see inset). The results correspond to 
Nei = 14 and = 27. 



energy level, separated by an energy gap from the rest of 
the spectrum (when compared with other energy separa- 
tions visible in the graph) . The gap becomes independent 
of the positions of the probes when the probe separation 
becomes larger than approximately 4 magnetic lengths. 
As the probes come closer to each other, the gap de- 
creases and, at the intersection points 6=0 and tt, the 
gap is completely closed, leading to an abrupt increase 
in the level degeneracy. 



V. BRAIDING THE MOORE-READ STATES 

A. General considerations. 

Let us consider the general situation when the number 
of flux has been increased by n units, N(f)=N^-\-n. We 
use the pinning Hamiltonian H{wi, . . . , W2n) [hq given in 
Eq. (18) to control the position of the anyons. Given 
that each Wj is a 2-dimensional variable, the Hamiltonian 



depends, parametrically, on 4n coordinates. We use the 
symbol x to denote these coordinates and the shorthand 
(x) for H{wi ,...,W2n)ino- 

Let us consider an arbitrary closed path in the parame- 
ter space: i x(-^) (xq =x(0)), where we use the length 
i to parametrize the loop. Now imagine that we move 
along this path with constant velocity v. One would like 
to study the time evolution U{t) generated by the time 
dependent Hamiltonian H{t) = H{-x{vt)): 



FIG. 6: (Color online) The figure illustrates a configuration 
in which all anyons are kept fixed except one, which is moved 
on the sphere. To compute different quantities at an arbitrary 
point P of the sphere, we introduce a local coordinate system 
{wi,W2) as described in the text. 



in the adiabatic limit, i.e. in the limit when the motion 

along the path is infinitely slow. 

If we denote by Px the projector onto the zero modes 

space of i?(x), then the classic Adiabatic Theorem says 
thatPtm 



U{t)P^,=Ua{t)P^,^0{v,t), 



(22) 



where Ua{t) is the adiabatic propagator, i.e. the unique 
solution of the fohowing system of equations;^ 



idtUait) = {H{t) + i[dtP^,P^])Ua{t), 

Ua{0) = I. 



(23) 



As long as the energy gap between the zero energy of 
the anyons and the energies of the excited states remains 
open, the adiabatic evolution is a good approximation 
of the true evolution, when restricted to the zero modes 
space. In fact, the two become identical if the velocity 
V is infinitely small. We want to point out that there 
are rigorous estimates (see for example Refs. [27H29|) of 
the errors o{v^t) that are made when one approximates 
U{t) by Ua{t)^ i.e. of the non-adibatic effects. Since the 
braids corresponding to the elementary TQC logic gates 
are quite long,'^ these estimates may play an important 
role in the field. 

The adiabatic evolution satisfies the following funda- 
mental relation: 



Px = f/a(x)PxoC^a(x)- 



(24) 



This equation tells us that the zero modes space is 
mapped into itself by the adiabatic evolution. This map- 
ping can be computed by solving Eq. (23), which tells 



idtU{t)=H{t)U{t), /7(0) = /, 



(21) 



us how the zero modes space is transported during the 
adiabatic evolution. One can put all these into a geo- 
metric perspective, by defining the following 1-form (the 



8 



adiabatic connection): 
1 



which is the same as Eq. (30). The numericahy calcu- 



da= -H{x.)d£^i[dP^,P^], 



(25) 



in which case Ua can be viewed as the solution of the 
fohowing differential equation: 



idUa = da Ua, t/a(xo) = /. 



(26) 



Thus, by fixing the speed we have eliminated the time. 
The quantum algorithms will use the adiabatic unitary 
transformations Ua{T), [/^(r'), . . . , resulted from taking 

the probes along certain closed paths F, F', 

Assume that, for each configuration of the anyons 
(equal to say that for each x), we chose a basis set (a 
gauge) V^i(x), ?/;/)(x) in the D=2^~^ dimensional zero 
mode space. We use '^(x) to denote the vector of compo- 
nents (x) , . . . , ^/^i) (x) . Due to the fundamental property 
of the adiabatic evolution Eq. (23), there exists a unitary 
DxD matrix W{x.) such that 

[/,(x)V^(xo) = iy(x)V^(x). (27) 

The Wilczeck-Zee connection is given byP^ 

di = iW'(x)-MW'(x), (28) 

which takes the classical form: 
1 



dA,,(x) 



-H{^)jid£-i{^j{^),dM^)), (29) 



where H{x.)ji are the matrix elements of the Hamiltonian 
in the chosen basis set. Since we are dealing with zero 
modes, the first term in the right hand side, above, is 
identically zero. 

The DxD unitary matrix W{x.) implements the adi- 
abatic evolution in the invariant subspaces Px'^o and it 
can be computed as the unique solution of the differential 
equation: 



zdW'(x) = W'(x) di,iy(xo) 



(30) 



which is to be integrated along the braiding path. For 
a loop F that starts and ends at xq, we solved this 
equation numerically, by considering a large number of 
points along the loop, xo,...,xx, and constructing the 
monodromy: 



Wr 



p p p p p 



(31) 



This amounts to finding the null space of the pin- 
ning Hamiltonian for each x/^. If we define Wr{k) = 
Pxo^x2---^xfc, then 



Wr{k + 1) - Wr{k) = Wr{k){P^,^, 

-Wr{k)P^,^,{P^,^, 

which is the finite difference version of 

idWr{x) = iWr{x)[dP^,P^]., 



Pyik)Pyik 

P^.), (32) 



lated Wr matrix becomes a unitary matrix only in the 
limit when the number of discrete points goes to infinity. 
To quantify how much does Wr deviate from a unitary 
matrix, we compute the absolute value of the determi- 
nant of Wr, which is compared to the value of 1, ap- 
propriate for a unitary matrix. In all the calculations 
presented in this paper, K was chosen large enough so 
that I det VFr 1=0.999 or better. This is a measure of how 
well converged are our numerical calculations. 



B. The Quantum Geometry of the zero modes 
states 



We can endow the zero modes with a curvature. The 
curvature form associated with the adiabatic connection 
is given by: 

dp = [d^A, - d,A^ - } dx^ A dx\ (34) 

The explicit expressions of its coefficients are: 

(F^,),,(x) = 2Im(S^V^,(x), [1 - Px]S.V^,(x)). (35) 

Besides the adiabatic connection and curvature, we can 
endow the parameter space with an intrinsic metric ten- 
sor, which we will refer to as the quantum metric ten- 
sor. First of all, we can introduce the following quantum 
distance!^ 



d^(x,xO = ||Px-Px'||//5, 



(36) 



where HS means the Hilbert-Schmidt norm. This dis- 
tance is at least second order differentiable in the coor- 
dinates X and x' and for this reason we can generate the 
quantum metric tensor ^^j, via a Taylor expansion: 



1 



d^(x,x + (5x) = -gl^{^)Sx^Sx^ ^o{S^'). 



(37) 



The coefficients of the quantum metric tensor are given 
by the classical expression!^ 

5«,(x) = 2Re^((9^V^(x), [1 - P^)]d,M^)). (38) 



In the following, we demonstrate an interesting relation 
between the adiabatic curvature and the quantum metric 
tensor. We fix all anyons, except one, in which case the 
parameter space becomes 2-dimensional. A point in this 
parameter space describes the position of the itinerant 
anyon on the sphere. To compute the coefficient of the 
curvature form at an arbitrarily chosen point P on the 
sphere, we introduce a local coordinate system by using 
the complex coordinate w = 2i?tan le"*"^, where (^, ^) 
are the usual spherical parameters when the North pole 
is fixed at P. We take wi=Re{w} and W2 =lm{w} as 
the two independent variables (see Fig. [6|. In these local 
coordinate system, the curvature form becomes 



(33) 



dF = F{w)dw^ A dw'^ 



(39) 
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FIG. 7: (Color online.) The adiabatic curvature F{P) ob- 
tained by moving one anyon while keeping the other fixed. 
Left-upper panel shows the results for odd number of elec- 
trons: N/N^=9/16, 11/20, 13/24, 15/28 and the right- 
upper panel shows the results for even number of electrons: 
Ar/Ar^=10/18, 12/22, 14/26, 16/30. The lower-left and lower- 
right panels show the sum and the difference between the 
odd and even results, respectively. For example, we added 
and subtracted the result for AT/AT^^IO/IS and N/N^=9/16, 
and then the results for N/N^=l2/22 and Ar/iV0=ll/2O, etc.. 




Area 

FIG. 8: (Color online.) The Berry phase accumulated by 
an anyon when moved along a path ^=const, with the other 
anyon fixed at the North pole. Different lines refer to differ- 
ent system sizes that are marked in the diagram with their 
corresponding N/N(f) numbers. A point on a particular line 
represents one braiding path, which is labeled by the area 
enclosed by the path, shown on the horizontal axis. 



where F is di DxD matrix. At P, in the same system of 
coordinates, we can prove the following general fact: 



gUP) = Tr{F(P)} g^P) 



(40) 



where is the standard metric of the sphere. This shows 
that the quantum metric and the standard metric of the 
sphere are related by a conformal transformation, which 
involves the trace of the curvature tensor. 

We start the proof of Eq. (40). The 2^~^ many body 



wavefunctions of the zero modes depend parametri- 
cally on the position w of the mobile anyon and take the 
following general form:^^ 



^^/i{w)=Ai{w)^^{w), z = l,...,2" 



(41) 



where the crucial fact is that the vectors ^i(^^;) are all 
analytic of w. Ai{w) are normalization factors that are 
not analytic of w. We have: 

Tr{F(^)} = 2Im^(Si^,(^), [1 - P^,]^2^^{w)). (42) 



Due to the presence of 1 — P^, we can ignore the action 
of the derivatives on the normalization constants Ai{w), 
which then can be pulled out of the scalar product as 
below: 



TrFM = 2El^iMI'x 

i 



(43) 



Next we use the fact that are analytic functions of 
i.e. they are functions of the single variable w rather 
than of two independent variables Wi and in which 
case: 



TTF{w) = 2Y^\Ai{w)\^x 

i 

lm{{d^^iiw), [1 - P^]d^^iiw))idiw)*d2w}. 



(44) 



The scalar product is real, which allows us to easily com- 
pute the imaginary part, resulting in: 



TvF{w) = 2Z\AiH\^x 

i 

(d^^iiw), [1 - P^]d^^i{w)). 



(45) 



If one repeats exactly the same arguments for the quan- 
tum metric tensor given in Eq. |38j the conclusion will be 
that: 



i 



(46) 



Since wi and W2 coincide with the geodesic coordinates of 
the sphere at P, i.e. the coordinates in which the metric 
tensor becomes the identity matrix when evaluated at P: 
g^{P)^jy = (5^^, the proof of Eq. 40 is completed. 
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VI. QUANTUM GEOMETRY FOR ONE FLUX 
ADDED. 

In this case we have 2 anyons and the zero modes space 
is 1-dimensional. We keep one anyon fixed and move the 
other along different braiding paths. The monodromy of 
any path T can be computed by integrating the curva- 
ture: 



jj^ i dF 



(47) 



where Sr is the surface enclosed by P. Thus, if we map 
the curvature, we can easily compute the monodromy of 
any arbitrary path. 

We compute the coefficient F of the curvature form 
(see Eq. 39) at a point P of the sphere from: 



Wr - I 

F{P) = hm 



(48) 



where F is a small path around P. The monodromy is 
computed as explained in the previous section and we 
obtain the limit by considering paths of decreasing ra- 
dius. The value of F computed this way coincides with 



the coefficient of the curvature in the local coordinate 
system {wi,W2) introduced in the previous section. The 
monodromy Eq. [48j however, is independent of the coor- 
dinates used to compute the curvature form. 

The upper panels in Fig. [t] plot F{P) as a function of 
the distance from P to the position of the fixed anyon. 
The two upper panels refer to odd and even numbers of 
electrons. The lower two panels show the sum and the 
difference between the results for odd and even number 
of electrons. Since one is interested in the thermody- 
namic limit, we plot sequences of curves for an increasing 
number of electrons. By comparing the curves for these 
sequences, one can determine how fast is the thermody- 
namic limit achieved. 

As one can clearly see, the curves in the top panels 
of Fig. [7| go asymptotically with the distance towards 
a constant value, which is precisely equal to the quasi- 
hole charge e*=e/4 (when working on the sphere, there 
is a small correction to this value, correction that goes 
to zero as the size of the sphere is increased). The most 
remarkable thing about Fig. [7] is that the curvatures for 
odd/even number of electrons have different thermody- 
namic limits, as one can clearly see by inspecting the dif- 
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FIG. 10: (Color online) Plots of the amplitude of the non-Abelina part of the curvature |f(P)| for different system sizes (each 
panel is marked with the corresponding N/Nc/)). For each size, we show |f(P)| calculated with the standard metric tensor (left) 
and with the quantum metric tensor (right). 



ference between the two, plotted in the lower-right panel. 

From the adiabatic curvature we compute the Berry 
phases accumulated (i.e. the exponent in Eq. 48) as we 



move one anyon along different paths of the form ^=ct., 
while keeping the other anyon fixed at the North Pole 
of the sphere. Fig. [8] shows the Berry phase as func- 
tion of the area enclosed by the path. We use this figure 
to draw several important conclusions. First, we point 
out that the Berry phase plotted in this figure includes 
also the Aharonov-Bohm phase due to the magnetic flux 
(I^B' i^AB=e'^<^B- Thus it is expected that the total Berry 
phase, as a function of the enclosed area, to go asymptot- 
ically to a linear curve, a feature that is obviously present 
in Fig. [Sj The slope of the asymptotic part of the curve 
is equal to the charge e*=e/4 of the anyons. 

Excepting the length of the lines, the graphs for dif- 
ferent system sizes look very similar, implying that the 
thermodynamic limit is achieved very fast. We expect the 
graphs for the larger systems to represent the thermody- 
namic limit with a high degree of accuracy. Looking at 
the largest systems, we can see that, for large enclosed ar- 
eas, the curves representing the Berry phase for systems 
with even number of electrons are shifted upward by tt 
relative to the curves representing the Berry phase for 
systems with odd number of electrons. This is a topo- 
logical effect, since the shift is independent of the area 
enclosed by the braiding loop or of the shape of the loop. 



as long as the loop is large enough. The finding is in 
total agreement with the SO spinorial representation of 
the braid group derived in Ref [7P| on the basis of under- 
ling CFT structure of the Moore-Read state. This effect 
was also considered the most direct signatu re of the Non- 
Abelian statistics for the Moore-Read state^^^ 

We conclude this subsection with the observation that 
the coefficient F depends on how we compute the area 



enclosed by F in Eq. 48, more precisely on what metric 
tensor is used when computing the area enclosed by the 
loops. The plots shown in Fig. [7| were obtained with the 
standard metric of the sphere. We have repeated the 
calculations using the quantum metric tensor Eq. (38) 
instead, in which case we have to re-scale the coefficient: 



F\w) = F{w)/^detg%{w). 



(49) 



According to the our previous analytic prediction (see 
Eq. 40), F^{P) should be identically 1. We numerically 



checked this prediction in the following way. The quan- 
tity that is readily available in the numerical calcualtions 
is the quantum distance. We can compute the determi- 
nant of the quantum metric tensor by considering a circle 
of radius p (in the standard metric of the sphere) centered 
at P and calculate the quantum distance between P and 
the points of this circle. If d\^ and d^ denote the max- 
imum, respectively minimum quantum distance to the 
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points of the circle, then 



detgl^{P) = hm 



(50) 



The computation of the determinant is done simultane- 
ously with the calculation of the curvature, which also 
requires the walk on the same circle. The numerical cal- 
culations give a direct confirmation that F^{P)=1 for all 
points of the sphere. 



VII. QUANTUM GEOMETRY FOR TWO 
FLUXES ADDED 

In this case we have 4 anyons and the zero modes space 
is 2-dimensional. We will keep 3 anyons fixed and move 
the forth one along different braiding paths. 

For t he No n-Abelian case, there is no simple Stokes 
theoremj ^^ l ^^ l which means the monodromy can not be 
simply computed from the curvature as we did for the 
previous case. Even so, mapping the curvature provides a 
clear picture of the non-comutative and topological prop- 
erties of the states. 

The parameter space remains 2-dimensional. As be- 
fore, a point in this parameter space indicates the posi- 
tion of the mobile anyon. Thus dF=Fdw^ Adw'^ ^ but F 
is now a 2x2 matrix. We compute F{P) using the same 
algorithm (see Eq. 48). Using the Pauli's matrices, a^, 
i=l, 2, 3, F{P) can be uniquely decomposed as: 



F{P) = fo{P) + f{P)-a, 



(51) 



where /o(P) = ^TrF(P) and f(P) is a 3-component vec- 
tor. We will refer to /o as the Abelian and to f-cr as 
the non- Abelian part of the curvature. It is important 
to notice how different quantities behave when changing 
the gauge, i.e. the basis in the 2-dimensional zero modes 
space. We have: fo{P) is gauge independent; the mag- 
nitude of f(P) is gauge independent; the orientation of 
f(P) is gauge dependent. 

Fig. [9] shows plots of TrP(P) for different system sizes. 
For each size, we show TrP(P) calculated with the stan- 
dard and with the quantum metric tensor. The numerics 
confirm again the theoretical prediction that TrP(P)=l. 

To demonstrate that the braid group is non- 
comutative, we need first to show that f(P) is non-zero. 
This, however, is not enough. We need also to rule out 
the existence of a particular gauge in which the adiabatic 
connection becomes diagonal at every point of the param- 
eter space. If such a gauge exists, the fiber bundle of the 
zero- modes degenerates into a trivial U{l)xU{l) fiber 
bundle, in which case there will be two 1-dimensional 
fibers that do not mix during the adiabatic braiding. 
Consequently, all monodromies will take a diagonal form 
in this gauge and they will commute with each other. 

As already mentioned, the magnitude of f(P) is gauge 
independent. Thus we can see if f(P) is zero or not by 




f(w+/Aw2) 



f(w+Awi) 



w+Awi 

FIG. 11: (Color online) A diagram of the parallel transport 
used in the calculation of the Twist density pTw- 



simply plotting its magnitude, which is shown in Fig. 10 
for different sizes. The graph clearly demonstrates that 
f(P) is non-zero and it appears to be concentrated near 
the positions of the fixed anyons. We will further discuss 



Fig. 10 in the next Section. 

Next, we introduce a scalar function which we call the 
Twist density pTw^ which gives a measure of how much 
are the fibers twisted during the adiabatic parallel trans- 
port. We start the construction from the following 2- 
form: 



dp = [D^F^ DvP] dw^ A dwjy^ 



(52) 



where I)^ denotes the covariant derivative correspond- 
ing to the adiabatic connection and [, ] denotes the usual 
commutator. This form is invariant to coordinate trans- 
formations, thus well defined. The coefficients of this 
form are 2x2 matrices. In the special case of a two 
dimensional parameter space, the form reduces to: 



dp = [DiF, D2F] dwi A dw2, 



(53) 



Based on the above observations, we construct the fol- 
lowing 2-form, 



dp^ 



det[P>iP,P>2P] dwi Adwi, 



(54) 



whose coefficient is a pseudo-scalar function. Since dp 
is invariant to coordinate transformations, dp^w is also 
invariant and hence well defined. The density 



PTw = ydet[DiF,D2F] 



(55) 



is gauge invariant and it is identically zero for trivial 
fiber bundles, in particular for U{l)xU{l) fiber bundle 
over our parameter space. Thus, if we show that Ptw is 
non-zero, that will be equivalent to demonstrating that 
the zero modes fiber bundle is non-trivial. 

Let us now give the physical interpretation of our con- 
struction. For this we consider three points on the sphere: 

w-\-Awi and w-\-iAw2^ as in Fig. 11 Assume that we 
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FIG. 12: (Color online) Plots of pTw{P) for different system sizes (each panel is marked with the corresponding N/N(j)). 
see how pTw{P) relates to the adiabatic curvature, for each size, we show |f(P)| on the left and pTw{P) on the right. 



To 



computed the curvature F at each of these points, in a 
pre-chosen, arbitrary gauge. The corresponding f (see 
Eq. 51) vectors are gauge dependent. To compare the 



three vectors, we need to adiabatically transport them to 
the same point w. For this we compute the monodromies 
Wi for the path w-\-Awi^w and W2 for the path 
w-\-iAw2^w and then compute Fi=WiF{w-\-Awi)Wi^ 
and F2=W2F{w-\-iAw2)W2'^- By decomposing Fi^2 as 
in Eq. 51, we obtain the parallel transport of f{w-\-Awi) 
and f(uH-iAw2) to w. We denote them by fi and f2, 
respectively. We now can ask if they are parallel. To 
quantify the answer to this question, we form the differ- 
ences ^if=fi-f and 62^=^2-^ and define 



Ptw = lim — 

Aw^O AWiAW2 



(56) 



This quantity is gauge invariant and measures how much 
are 5if and 82^ deviating from being parallel. The two 
definitions of p^w given in Eqs. (55) and (56) coincide. 



Fig. 12 shows plots of the Twist density for different sys- 
tem sizes. 



VIII. DISCUSSION 

We now can draw several conclusions regarding 
the non-commutative and topological properties of the 
Moore-Read states with four anyons. Although the size 



of the systems we attempted are still small, the trend seen 
in the sequence of plots shown in Fig. [lO] suggests that 
the Non-Abelian curvature is localized in the vicinity of 
the fixed anyons. This implies a topological property for 
the Non-Abelian part of the mononodromies, i.e. their 
independence of the shape of braiding paths as long as 
the braiding path are far enough from the anyons. 



Fig. 12 demonstrates the existence of non-commutative 
monodromies. Indeed, if all the monodromies were com- 
muting with each other, that will imply the existence 
of a gauge in which the fiber-bundle of the zero modes 
becomes trivial. But this will imply Ptw=0, which is 



contradicted by Fig. 12 



Fig. 12 reveals much more, namely the splitting of 
the zero modes space as predicted by the fusion rules 
of the underlying CFT structure of the Moore-Read 
sequence. I^^Q^ For this, notice that the Twist density is 
practically zero when the mobile anyon comes near the 
fixed anyons and takes non-zero values only in between 
the fixed anyons. This implies that, although |f | takes ap- 
preciable values in the regions near the fixed anyons, the 
states are commutative. In other words, when the mobile 
anyon comes close to a fixed anyon, the two eigenmodes 
of the curvature matrix F{P) split the zero modes space 
into two sectors that don't mix with each other during 
the braiding. If things happen as predicted by the fusion 
rules of the underlying CFT structure of the Moore- Read 
sequence, then we should be able to read from the classic 
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FIG. 13: (Color online.) Comparison between the results 
obtained for the Abelian case (one flux added) and for the 
non-Abelian case (two flux added). The left/right panels 
refer to odd/even number of electrons. The solid lines in 
the upper/lower panels is the difference/sum between the 
Abehan Berry curvature for Ar/iV^=16/30 and Ar/iV0=15/29 
see Fig.|9|. The two open circle lines in the upper/lower pan- 
els show the values (multiplied by 2) of /o(x)/|f(x)| recorded 
along the path shown in the inset. In the left panels, the dot- 
ted lines correspond to the system sizes A/'/A^0= 15/29 and 
13/25. In the left panels, the dotted lines correspond to 
iV/Ar^=14/27 and 12/23. 



Bratelli diagram what these sectors are: they should be 
the q-spin S'=0 and 5'=1 sectors.!^ For the sphere geom- 
etry, these sectors correspond to the even, respectively 
the odd numbers of electrons. 

To demonstrate that this is the case, we notice that 



|f(P)| plotted in Fig. 10 is equal to half the difference 
between the two eigenvalues of the curvature coefficient 
F{P) and that TrF(P) plotted in Fig. [g] is equal to the 
sum of the two eigenvalues. Thus, near the fixed quasi- 
holes, the values of these quantities, properly normalized, 
should compare well with the values shown in the plots of 
the lower panels of Fig. jT] To verify that this is the case, 
we generated line-plots out the surface plots shown in 
Figs. [10] and [9| The line plots were generated by record- 
ing the points along the path shown in the inset of Fig. 13 
We used the data corresponding to the standard metric 
of the sphere. These line plots are compared in Fig. [13] 
with the absolute value of the data shown in Fig. [7| As 
one can see, the matching is almost perfect near the fixed 
anyon. By comparing two system sizes, we can also see 
the agreement becoming better as the size of the system 
is increased. We expect the agreement to become perfect 
as the sphere radius is taken to infinity. It is important 
to notice that there is agreement for both even and odd 
number of electrons cases. 

Once we made this connection, we can we can go back 



and discuss the localization of the non-abelian curva- 
ture near the fixed quasiholes. The difference between 
the Berry phases for odd and even number of electrons, 
shown in Fig. [sj converges extremely fast to tt [the plot 
suggests an exponential convergence]. From the connec- 
tion made in the previous paragraph, this implies that 
the non-abelian curvature is also strongly [exponentially] 
localized near the fixed anyons. Unfortunately, the size 
of our systems is too small to allow a quantitative eval- 



uation of this asymptotic decay behavior. Fig. [13] tells 
us how far are we from a converged situation, where the 
open circles should overlap with the continuous line. 

We end this discussion Section by mentioning that we 
did compute several braiding monodromies and tried to 
compare their group properties under multiplication with 
the predictions following from the underlying CFT struc- 
ture of the Moore- Read state derived in Ref. [191 Unfor- 
tunately, the size of our system is too small to see a 
strong correlation between the two. This can be seen 
from Fig. [13] where the continuous line practically repre- 
sent the thermodynamic limit of the non-abelian curva- 
ture and the open circles represent the same quantity for 
our largest finite system. The most problematic part is 
the large discrepancy seen at points in between the quasi- 
holes (the end of the curves in Fig. 13), which affects 



any monodromy looping around a fixed anyon. 



IX. CONCLUSIONS 

We have shown that the Hall sequences can be viewed 
as the zero modes of certain Hamiltonians, which were 
written down explicitly using many-particle creation op- 
erators. The compressible Hall states with the anyons 
fixed at definite locations were shown to be the zero 
modes of a pinning Hamiltonian, which was also written 
down explicitly. We have developed an efficient diago- 
nalization algorithm for the pinning potential, which was 
subsequently used to explore some general properties and 
to map the quantum geometry of the Moore-Read states 
with two and four anyons. 

Working with two anyons, we were able to give the 
first direct confirmation of the topological properties of 
the monodromies. The monodromy corresponding to a 
loop enclosing one of the anyons was found to differ by 
exactly a factor -1, when computed for even/odd num- 
ber of electrons (which correspond to the q-spin 5'=0/l 
sectors). This is the first explicit confirmation of the the- 
oretical predictions based on the underlying CFT struc- 
ture of the Moore-Read sequence. 

By fusing the two anyons, we found that the electron 
density is different for the even/odd number of electrons. 
This is the first direct confirmation that one can deter- 
mine if a quantum state is in the 5=0/1 sectors by fusing 
the anyons and measuring the electron density. This type 
of measurement stands at the basis of read-in and read- 
out processes of TQC.^^ 

Working with four anyons, we mapped the abelian and 
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non-Abelian parts of the adiabatic curvature for the case 
of one itinerant anyon and three anyons fixed at the equa- 
tor, in the most possible spread configuration. The map- 
ping reveals that the non-Abelian curvature is strongly 
localized near the fixed quasiholes. We introduced the 
Twist density, which measures twisting of the zero modes 
during the adiabatic braiding. We found that the Twist 
density is practically zero near the fixed anyons, fact that 
signaled a splitting of the zero modes in two non-mixing 
sectors. Further analysis showed that this splitting is 
precisely the one implied by the fussion rules of the un- 
derlying CFT structure of the Moore-Read sequence. 
If the computations can be implemented to larger sys- 



tems sizes, the present study open the possibility of: sim- 
ilar studies for higher level Hall sequences that support 
universal quantum computation; direct implementation 
and verification of the quantum gates found in Ref. |T71 
and a direct simulation of a quantum algorithm, includ- 
ing the read in, braiding and read out phases. 
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